PL.TR-93-2192 


AD-A278  137 

■llllllll 


INVESTIGATION  OF  LATERAL  VARIATION  IN  THE  SEISMIC 
VELOCITY  STRUCTURE  OF  THE  SHALLOW  CRUST  BENEATH 
EASTERN  MASSACHUSETTS  AND  SOUTHERN  NEW  HAMPSHIRE 


Susan  E.  D'Annolfo 
Alan  L.  Kafka 


Weston  Observatory 

Department  of  Geology  and  Geophysics 

Boston  College 

Weston,  MA  02193 


30  September  1993 


Scientific  Report  No.  2 


Approved  for  public  release;  distribution  unlimited 


94-11550 

J I  mill  mil  mil  mil  mil  mil  li  llll  lill 


PHILLIPS  LABORATORY 

Directorate  of  Geophysics 

AIR  FORCE  MATERIEL  COMMAND 

HANSCOM  AIR  FORCE  BASE,  MA  01731-3010 


The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  representing  the 
official  policies,  either  expressed  or  implied,  of  the  Air  Force 
or  the  U.S.  Government. 

This  technical  report  has  been  reviewed  and  is  approved  for 
publication. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OM8  No  0704  0188 

*0'  t*i*i  <o<'en>oA  of  infortw^tioi*  to  *  *»Ow/  •nciudm^  t^e  fo'  ifntru<tio^.  ea>tiin9  o«t«  >ou'cei. 

«nO  m«tr«t««Ainq  0«t«  n««40<l.  *^0  (0<*'0(«t«A9  Sf*0  the  <Oltf<t»OA  O*  iMor>n*tion  S^nd  COmmonti  rMA'Oin^  thty  burden  etttm«te  or  eny  ot^er  «»oe<t  of  tn<« 

<0«ect‘0n  0*  >nto*n«ation  -nOwtfiA^  tw99ett'Ony  *0*  'ediKtn9  in«%  burden  to  'AreitifAqion  HeedQueriers  Serwi<e).  0«rectOrete  'o*  >ntorm«tiOn  Operation^  and  Aeportt.  w  IS  Jener>on 
DavtM.fhMcy.Su'te  i?04  anm^ton.  va  2^^oiA3Q2.  and  to  the  OTfue  o*  Mane9entent  and  •wd9et  .  Paperwork  KedvKttonerotect(0r04^)tM}.  vva%n<A9ton.  OC  20S03 

1.  AGENCY  USE  ONLY  (Ltivt  blink)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

30  September  1993  Scientific  No.  2 

4.  TITLE  AND  SUBTITLE 

Investigation  of  Lateral  Variation  in  the  Seismic  Velocity  Structure  of  the 
Shallow  Crust  Beneath  Eastern  Massachusetts  and  Southern  New  Hampshire 

S.  FUNDING  NUMBERS 

PE  62101F 

PR  7600  TA  09  WU  AL 

Contract  F19628-90-K-0035 

6.  AUTHOR(S) 

S.E.  D'Annolfo  and  A.L.  Kafka 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORE  SS(ES) 

Weston  Observatory 

Department  of  Geology  and  Geophysics 

Boston  College 

Weston,  MA  02193 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING /MONITORING  AGENCY  NAM£(S)  AND  AODRESS(ES) 

Phillips  Laboratory 

29  Randolph  Rd. 

HanscomAFB.MA  01731-3010 

Contract  Manager  James  Lewkowicz/GPEH 

10.  SPONSORING /MONITORING 

AGENCY  report  NUMBER 

PL-TR-93-2192 

11.  SUPPLEMENTARY  NOTES 

1Z«.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

IZb.  DISTRIBUTION  CODE 

13  ABSTRACT  (MiMimum  200  wordi) 

This  is  one  of  five  scientific  reports  describing  spiecific  research  projects  conducted  at  Weston 
Observatory  under  Contract  No.  Fl%28-90-K-0035.  The  research  conducted  under  this  contract 
covers  a  range  of  topics  related  to  seismology  in  general  and  to  nuclear  test  monitoring  in 
particular.  'Oiis  report  describes  a  study  in  which  group  velocity  dispersion  was  determined  for 
short-period  Rayleigh  waves  recorded  from  blasts  detonated  at  the  San-Vel  quarry  in  Littleton, 

MA.  Field  data  were  recorded  to  complement  network  data  in  order  to  create  a  denser  distribution 
of  stations  than  was  possible  in  past  experiments.  Group  velocity  dispersion  was  estimated  for 
each  path.  The  resulting  dispersion  curves  were  analyzed  to  investigate  the  relationship  between 
group  velocities  and  sh^low  crustal  stracture  beneath  eastern  MA  and  southern  NH.  Our  results 
suggest  that  the  quarry  lies  in  a  "trough"  of  low  group  velocities,  where  the  group  velocities 
systematically  increase  toward  the  E  and  W-NW  directions.  This  creates  a  "valley"  of  group 
velocities  where  the  primary  trend  of  the  valley  trough  is  N,  and  a  secondary  feature  is  a  trough 
that  trends  in  a  NE  direction.  This  NE  trending  feature  is  of  particular  interest  because  of  a 
possible  correlation  with  the  trend  of  the  Clinton-Newbury  fault  zone. 

14  SUBlEa  TERMS 

Shallow  crustal  structure,  Rg  wave  dispersion,  quarry  blasts, 
chemical  explosions 

IS.  NUMBER  OF  PAGES 

104 

16.  PRICE  CODE 

17.  security  classification  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  Unclassified 

ZO.  LIMITATION  OF  ABSTRACT 

SAR 

MSN  7S40-01-280-SS00  Standard  Form  298  (Rev  2-89) 


by  AN^  Std 


11 


CONTENTS 


1.  INTRODUCTION 

t .  1  Background  and  Purpose  of  This  Study  1 

1.2  Data  Acquisition 

1.3  Seismic  Sources 

1.4  Instrumentation 

1.5  Description  of  Data  Set  1 

2.  GEOLOGY  AND  CRUSTAL  STRUCTURE  OF  SOUTHERN  NEW  ENGLAND 

2.1  Geology  of  Southern  New  England  21 

2.2  Crustal  Structure  of  Southern  New  England  24 

3.  SURFACE  WAVE  THEORY 

3. 1  Rg  Waves  28 

3.2  Measuring  Group  Velocities  29 

4.  DATA  PROCESSING 

4.1  Origin  Times  31 

4.2  Effects  of  Origin  Time  Errors 

on  Group  Velocity  Measurements  36 

5.  RESULTS 

5.1  Summary  of  Results  From  Previous  Studies  of  SNE  40 

5.2  Results  of  This  Study  42 

6.  DISCUSSION  AND  CONCLUSIONS  61 


REFERENCES 
Appendix  A: 
Appendix  B: 


Group  Velocity  Data 
Group  Velocity  Curves 


62 

65 

75 


iii 


^  «o  «o  a> 


IV 


PREFACE 


This  is  one  of  five  scientific  reports  describing  specific  research  projects  conducted  at  Weston 
Observatory  under  Contract  No.  F19628-90-K-0035.  The  research  conducted  under  this  contract 
covers  a  range  of  topics  related  to  seismology  in  general  and  to  nuclear  test  monitoring  in 
particular. 

This  scientific  report  consists  of  an  M.S.  thesis  written  by  Susan  E.  D'Annolfo  under  the 
supervision  of  Professor  Alan  L.  Kafka.  In  this  study,  group  velocity  dispersion  was  determined 
for  Rayleigh  waves  between  periods  of  0.2  and  2.2  sec  (Rg)  recorded  from  blasts  detonated  at  the 
San-Vel  quarry  in  Littleton,  MA.  Field  data  were  recorded  to  complement  netwoiit  data  in  order  to 
create  a  denser  distribution  of  stations  than  was  possible  in  past  experiments.  Rg  wave  dispersion 
was  estimated  for  each  path.  The  resulting  dispersion  curves  were  analyzed  to  investigate  the 
relationship  between  Rg  wave  velocities  and  the  shallow  crustal  structure  beneath  eastern 
Massachusetts  and  southern  New  Hampshire. 

The  results  of  this  study  suggest  that  there  is  a  systematic  lateral  variation  in  the  structure  of 
the  shallow  crust  underlying  eastern  Massachusetts  and  southern  New  Hampshire.  Lateral 
variation  in  group  velocities  widiin  this  area  appear  to  depend  on  distance  from  the  San-Vel  quarry. 
Three  dimensional  plots  of  the  group  velocity  data  reveal  that  the  San-Vel  quarry  appears  to  lie  in  a 
"trough"  of  particularly  low  group  velocities,  where  the  group  velocities  systematically  increase 
toward  the  east  and  west-northwest  directions.  This  creates  a  "U-shaped  valley"  of  group 
velocities  where  the  primary  trend  of  the  valley  trough  trends  in  a  north-south  direction,  and  a 
secondary  feature  is  a  "valley"  that  trends  in  a  northeast  direction.  This  northeast  trending  featureis 
of  particular  interest  because  of  a  possible  correlation  with  the  structural  geology  of  the  area, 
particulary  the  trend  of  the  Clinton-Newbury  fault  zone,  which  forms  the  boundary  between  the 
Merrimack  trough  and  the  Putnam-Nashoba  terrane. 


V 


INVESTIGATION  OF  LATERAL  VARIATION  IN  THE 
SEISMIC  VELOCITY  STRUCTURE  OF  THE  SHALLOW  CRUST  BENEATH  EASTERN 
MASSACHUSETTS  AND  SOUTHERN  NEW  HAMPSHIRE 


1.  INTRODUCTION 

1.1  Background  and  Purpose  of  this  Study 

Studies  of  shallow  crustal  structure  using  Rg  wave  dispersion  have  shown  that  the 
dispersion  patterns  in  Southern  New  England  (SNE)  exhibit  differences  in  group  velocity  from 
one  area  to  the  next  (e.g.  Kafka  and  Doilin,  1985;  Kafka,  1988).  Rg  is  a  short  period 
fundamental  mode  Rayleigh  wave,  and  the  dispersive  properties  of  Rg  waves  are  sensitive  to 
variations  in  the  seismic  velocity  structure  of  the  shallow  crust.  Based  on  observed  differences 
in  Rg  dispersion,  Kafka  and  Skehan  (1990)  divided  SNE  into  five  regions  of  distinct  Rg 
dispersion  characteristics  that  they  called  "dispersion  regions."  The  most  clearly  distinct 
dispersion  region  in  SNE  is  the  northern  portion  of  the  Hartford  Rift  Basin  (Rgure  1),  where 
Rg  group  velocities  are  distinctly  Tower  than  they  are  in  other  parts  of  SNE  (presumably  due  to 
the  abundance  of  sediments  and  sedimentary  rock  in  that  area).  On  the  other  hand,  systematic 
lateral  variations  of  Rg  group  velocities  have  not  been  clearly  delineated  in  crystalline 
basement  provinces.  One  crystalline  basement  province  that  was  thought  to  be  characterized  by 
distinctly  higher  Rg  group  velocities  was  southwestern  Connecticut  (Kafka  &  Doilin,  1985).  In 
a  more  recent  study  however,  Kafka  and  Bowers  (1991)  have  shown  that  southwestern 
Connecticut  is  actually  characterized  by  Rg  group  velocities  quite  similar  to  those  of  other  areas 
where  the  crystalline  basement  rocks  are  at  or  near  the  surface.  Thus  the  question  remains: 
"What  is  the  pattern  of  lateral  variation  of  shallow  crustal  structure  in  the  crystalline 
basement  rock  underlying  SNE?" 
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Figure  1 :  Map  of  Rg  dispersion  regions  in  Southern  New 
England,  (after  Kafka  and  Skehan,  1990) 


The  Rg  dispersion  pattern  in  the  rather  large  area  that  extends  from  the  Bronson  Hill 
Anticlinorium  to  the  Avalonian  Terrane  (Figure  2)  displays  very  little  variation  in  group 
velocities.  This  observation  suggests  that  the  shallow  crust  underlying  this  area  is  laterally 
homogeneous  (at  least  at  the  scale  of  features  revealed  by  previous  Rg  dispersion  studies). 

Because  of  that  lack  of  observed  lateral  variation,  Kafka  and  Skehan  (1990) 
characterized  that  entire  area  as  an  Rg  dispersion  region,  which  they  called  the  Bronson  Avalon 
Dispersion  Region  (BADR).  Since  the  BADR  includes  a  wide  range  of  geoiogicat  structures  and 
rock  types,  it  is  surprising  that  this  region  appears  to  be  so  homogeneous.  The  Rg  dispersion 
pattern  seems  to  ignore  the  transition  from  one  geological  feature  to  the  next.  In  contrast,  there 
is  at  least  one  other  area  of  New  England  (southeastern  Maine)  where  significant  differences  in 
Rg  dispersion  have  been  observed  that  appear  to  correlate  with  the  geology  of  the  crystalline 
basement  structures.  Kafka  and  Reiter  (1987)  found  evidence  for  lateral  anisotropy  in  the 
shallow  crust  of  southeastern  Maine,  where  the  trend  of  the  anisotropy  is  parallel  to  the 
structural  grain  of  the  Appalachians.  Thus  the  question  arises:  Does  the  seismic  velocity 
structure  of  the  shallow  crust  beneath  the  so  called  BADR  really  have  no  systematic  relationship 
to  the  surface  mapped  geology?  AHematively,  does  the  relationship  between  the  surface  geology 
and  shallow  crustal  structure  exist  at  a  scale  that  is  too  small  to  be  seen  by  the  distribution  of 
paths  described  in  the  above  mentioned  studies? 

During  the  summer  of  1987  an  experiment  was  performed  where  portable  seismic 
stations  were  set  up  around  the  San-Vel/Lonestar  Quarry  in  Littleton,  MA.  Since  the  stations 
were  densely  distributed  and  the  distances  between  source  and  receiver  varied  from  a  few 
kilometers  to  71  km  (Figure  3).  this  experiment  provided  an  opportunity  to  investigate  the 
pattern  of  Rg  dispersion  within  the  BADR  in  greater  detail  than  previously  possible.  Figure  4 
shows  group  velocities  measured  from  six  seismograms  recorded  during  the  1987  experiment. 
In  spite  of  this  denser  station  coverage,  there  is  still  very  little  variation  in  the  group 
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Figure  2:  Map  of  tectonic  regions  in  Southern  New  England,  with  the 
BAOR  shaded,  (after  Kafka  and  Skehan,  1990) 
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velocities.  This  is  surprising  considering  that  the  propagation  paths  traverse  geotogicaliy 
diverse  areas. 

It  is  possible  that  systematic  velocity  variations  actually  do  exist  within  the  BAOR  but 
are  not  evident  because  the  seismic  waves  that  have  been  obsenred  there  have  propagated  aaoss 
relatively  large  (^stances.  Source-to-receiver  distances  in  previous  studies  have  ranged  from  a 
little  less  than  25  km  to  over  100  km,  and  when  an  Rg  wave  travels  a  great  (fistance  it  travels 
with  a  group  velocity  that  reflects  an  average  of  the  different  types  of  materials  it  encounters. 

For  my  thesis  research,  I  investigated  the  details  of  lateral  variations  in  seismic 
velocity  structure  of  the  shallow  crust  beneatii  eastern  MA  and  southern  NH  in  greater  detail 
than  has  been  done  in  the  past  by  utilizing  field  data  to  obtain  a  denser  station  spacing.  The 
object  of  this  experiment  is  twofold:  1)  to  obtain  paths  of  shorter  iength;  and  2)  to  have  a  more 
dense  distribution  of  paths.  The  purpose  of  this  investigation  is  to  determine  whether  or  not  the 
lack  of  obsen/ed  variation  of  Rg  dispersion  in  SNE  is  an  artitect  of  the  relatively  sparse 
distribution  of  paths. 

In  this  study,  the  following  two  alternative  hypotheses  regarding  the  extent  of  lateral 
variations  in  seismic  velocity  structure  of  the  shallow  crust  beneath  eastern  MA  and  southern 
NH  were  tested: 

(1)  The  crystalline  basement  rock  underlying  eastern  MA  and  southern  NH  actually  is 
homogeneous  at  the  scale  of  features  differentiated  by  Rg  waves,  even  with  the  addition  of  the 
field  stations. 

(2)  Lateral  variations  exist  within  the  shallow  crust  underlying  eastern  MA  and  southern  NH, 
but  they  exist  on  a  scale  that  is  too  small  to  have  been  revealed  by  previous  studies. 

As  will  be  discussed  below,  the  results  of  this  study  suggest  that  there  does  in  fact  appear 
to  be  smaller  scale  lateral  variation  in  the  shallow  crust  beneath  the  area  surrounding  the  San- 
Vel  quarry. 
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NA  =  Proto-North  Amen  can  Terrene 
CV  =  Connecticut  Valley  Syncllnorlum 
CL  =  Cameron's  Line 
HRB  =  Hartford  Rift  Basin 
BH  =  Bronson  Hill  AntlcHnoiium 


KC  =  Kearsarge-Central  ME 
Syncllnorlum 
MT  =  Merrimack  Trough 
PN  =  Putnam-Nashoba  Terrene 
AV  =  Avalonian  Superterrane 


Figure  3:  Tectonic  map  of  Southern  New  England  showing  the  sites  of 
the  1987  Littleton  Quarry  Blast  Experiment. 
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1 .2  Data  Acquisition 

For  this  research,  analysis  was  done  on  the  group  velocity  dispersion  of  Rg  waves 
recorded  from  quarry  blasts  detonated  at  the  San-Vel/Lonestar  Quarry  in  Littleton,  MA  in  order 
to  investigate  the  extent  of  lateral  variation  in  the  shallow  crust  beneath  the  area  surrounding 
that  quarry.  There  are  many  seismograms  of  quarry  blasts  and  shallow  focus  earthquakes  in  the 
New  England  area  that  have  been  recorded  digitally  by  the  New  England  Seismic  Network  (NESN) 
operated  by  Weston  Observatory.  The  Earth  Resources  Laboratory  at  the  Massachusetts 
institute  of  Technology  (MIT)  also  operates  a  regional  network  that  records  seismic  events  in 
this  region.  Data  from  both  of  these  networks  were  used  for  this  study.  Both  networks  record 
clear  signals  and  prominent  Rg  waveforms  from  the  San-Vel  quarry.  Data  from  several  field 
experiments  performed  by  the  Air  Force  Geqshysics  Laboratory  (AFGL)  at  Hanscom  AFB,  that 
were  part  of  the  1987  field  study  described  above  are  also  included  in  this  study.  In  addition, 
portable  seismic  stations  were  installed  in  the  field  to  investigate  whether  or  not  a  denser 
station  spacing  would  reveal  any  lateral  variations  in  seismic  velocity  structure  that  have  yet 
been  observed  in  previous  studies. 

As  part  of  this  study,  a  field  experiment  was  designed  in  which  portable  field  stations 
were  set  up  to  follow  a  straight  line  path  that  extended  from  the  San-Vel  Quarry  out  in  a 
northwest  direction  to  New  Ipswich,  NH.  The  path  was  approximately  35  km  long.  Eleven 
receiver  stations  were  set  up.  The  ot^ect  was  to  locate  the  stations  at  intervals  ranging  from  2 
km  to  4  km  along  the  line.  Most  of  the  stations  were  located  on  private  property  where 
permission  was  granted  from  the  owners.  A  few  of  the  sites  fell  on  state  or  city  owned  property. 
Permission  was  also  obtained  from  the  proper  authorities  to  utilize  those  areas.  Sites  were 
numbered  in  a  series  labeled  NW01  through  NW1 1  to  delineate  the  line  direction  and  station 
number.  Of  the  eleven  original  field  sites,  only  one  (NW09)  was  abandoned  because  it  was  too 
difficult  to  access  on  short  notice.  Equipment  was  deployed  at  the  remaining  ten  sites  once  per 
week  during  the  months  of  July  and  the  beginning  of  August,  1989.  Because  of  technical 
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failures  and  human  error,  data  recovery  for  the  experiment  was  approximately  65%  successful 
for  the  five  blasts  undertaken  to  record. 

There  was  an  attempt  to  include  in  this  study  a  second  and  third  line  of  field  stations  that 
extended  out  in  a  northeast  and  a  southeast  direction.  However,  time  only  allowed  me  to  record  a 
small  part  of  the  northeast  line  once  before  all  of  the  field  equipment  became  unavailable  for  the 
remainder  of  the  season.  Figure  5  is  a  map  showing  ail  of  the  paths  used  for  this  study  and  their 
proximity  to  the  San-Vel  Quarry  (located  at  42.554'*N  and  71.517**W).  Figure  6  is  a  diagram 
of  the  quarry  pit  which  encompasses  approximately  one  square  kilometer  of  land.  Because 
different  faces  of  the  pit  are  mined  at  cfifferent  times,  an  average  location  was  chosen  at  the 
center  of  the  quarry  pit  to  represent  the  geographic  location  of  the  quarry.  This  means  that  the 
location  of  any  particular  blast  could  at  worst  be  off  by  one  half  of  a  kilometer.  Table  1  lists  all 
of  the  stations  by  name  with  their  geographical  coordinates  and  distances  from  the  quarry. 

1.3  Seismic  Sources 

Quarry  blasts  detonated  at  the  San-Vel  Quarry  were  chosen  to  be  the  primary  seismic 
sources  for  this  research  for  several  reasons:  1)  the  quarry  operator  was  extremely 
cooperative  and  worked  cioseiy  with  me  on  the  experiment,  enabling  me  to  pinpoint  source 
location  and  quite  accurately  determine  event  origin  times;  2)  blasts  from  the  quarry  are  often 
large  and  generate  prominent  Rg  waveforms;  3)  there  are  two  to  three  quarry  blasts  per  week 
at  this  quarry  during  summer  months;  and  4)  there  was  a  source  of  AFGL  data  from  this  quarry 
from  the  1987  experiment  that  augmented  fois  data  set  quite  well  and  had  not  been  analyzed  yet. 

1 .4  Instrumentation 

The  NESN  stations  operated  by  Weston  Observatory  consist  of  1  Hz  HS-10  vertical 
component  sensors,  and  the  MIT  network  stations  consist  of  1  Hz  Mark  Products  L-4C  vertical 
component  sensors.  Data  are  telemetered  to  Weston  and  MIT  by  telephone  line  and  are  recorded 
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Figure  6:  Map  of  5an-Ve1  Quarry  mining  pit.  The  star  in  the  center  of  the 
map  delineates  the  geographic  location  of  the  quarry  used  for  data 
processing  as  described  in  the  text. 


Table  1 


Station 

Lai 

Lon 

NWOl 

42.5790 

-71.5800 

NW02 

42.5940 

-71.5960 

NW03 

42.5950 

-71.6110 

NW04 

42.6155 

-71.6630 

NH05 

42.6240 

-71.6830 

NW06 

42.6360 

-71.7010 

NW07 

42.6520 

-71.7250 

NW08 

42.6610 

-71.7540 

NWIO 

42.6880 

-71.8880 

NWll 

42.7550 

-71.8500 

WES 

42.3847 

-71.3221 

QUA 

42.4566 

-72.3738 

WFM 

42.6106 

-71.4906 

ONH 

43.2792 

-71.5056 

PNH 

43.0942 

-72.1358 

WKH 

43.8683 

-71.3997 

DNH 

43.1225 

-70.8948 

GLO 

42.6403 

-70.7272 

UXB 

42.0614 

-71.6773 

DUX 

42.0686 

-70.7678 

COD 

41.6858 

-70.1350 

AZl 

42.3864 

-71.5311 

AZ2 

42.3937 

-71.6467 

AZ3 

42.5156 

-71.8221 

AZ4 

42.6079 

-71.7821 

LCl 

42.5312 

-71.5571 

LC2 

42.5071 

-71.5910 

LC3 

42.4962 

-71.6267 

LC4 

42.4603 

-71.6935 

P102 

42.5400 

-71.5660 

P103 

42.5360 

-71.5880 

P107 

42.5300 

-71.6900 

PllO 

42.5200 

-71.8000 

P112 

42.5145 

-71.8396 

P113 

42.5121 

-71.8546 

P115 

42.5187 

-71.9054 

P201 

42.5480 

-71.5390 

P202 

42.5450 

-71.5460 

P203 

42.5440 

-71.5550 

P206 

42.5310 

-71.5980 

NSAO 

42.4400 

-71.7270 

NSAl 

42.4390 

-71.7280 

NSA2 

42.4400 

-71.7240 

NSA3 

42.4410 

-71.7280 

NSBl 

42.4400 

-71.7290 

NSB2 

42.4420 

-71.7260 

NSB3 

42.4390 

-71.7230 

NSB4 

42.4370 

-71.7260 

NSB5 

42.4390 

-71.7300 

P»s«  Affiliation 


5.9 

FIELD  1989 

7.9 

FIELD  1989 

9.0 

FIELD  1989 

13.8 

FIELD  1989 

15.7 

FIELD  1989 

17.6 

FIELD  1989 

20.3 

FIELD  1989 

22.8 

FIELD  1989 

33.9 

FIELD  1989 

35.3 

FIELD  1989 

25.0 

BC  NETWORK 

71.0 

BC  NETWORK 

6.7 

MIT  NETWORK 

80.6 

MIT  NETWORK 

78.5 

MIT  NETWORK 

146.3 

MIT  NETWORK 

81.1 

MIT  NETWORK 

65.5 

MIT  NETWORK 

56.3 

MIT  NETWORK 

82.0 

MIT  NETWORK 

149.5 

MIT  NETWORK 

18.4 

BC/MIT  1987 

20.4 

BC/MIT  1987 

25.1 

BC/MIT  1987 

22.4 

BC/MIT  1987 

3.8 

BC/MIT  1987 

7.7 

BC/MIT  1987 

10.7 

BC/MIT  1987 

17.5 

BC/MIT  1987 

4.0 

AFGL  1987,88 

5.9 

AFGL  1987,88 

14.0 

AFGL  1987 

23.1 

AFGL  1987 

26.9 

AFGL  1987 

27.8 

AFGL  1987 

31.9 

AFGL  1987 

1.5 

AFGL  1988 

2.1 

AFGL  1988 

2.9 

AFGL  1988 

6.6 

AFGL  1988 

21.5 

AFWL  1987 

21.6 

AFWL  1987 

21.3 

AFWL  1987 

21.5 

AFWL  1987 

21.6 

AFGL  1987 

21.2 

AFGL  1987 

21.2 

AFGL  1987 

21.6 

AFGL  1987 

21.7 

AFGL  1987 
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digitally  at  50  samples  per  second  (sps)  at  Weston  and  100  sps  at  MIT,  when  a  trigger  oorufition 
is  met. 

Equipment  for  most  of  the  field  stations  was  supplied  by  the  Air  Force  Geophysics 
Laboratory  (AFGL)  at  Hanscom  AFB,  Bedford,  MA.  Each  of  these  stations  consisted  of  a  2  Hz, 
three  component  Sprengthenether  S-6000  seismometer  connected  to  a  Terra-Technology  DCS- 
302  digital  recorder.  Some  of  the  field  stations  consisted  of  portable  instruments  from  Weston 
Obsenratory.  These  stations  consisted  of  a  2  Hz,  three  component  Sprengthenether  S-6000 
seismometer  connected  to  a  Sprengthenether  DR-200  digital  recorder.  All  field  recorder  docks 
were  synchronized  with  a  Kinemetrics  satellite-synchronized  clock  both  before  deployment  and 
upon  return  from  the  field  in  order  to  compensate  for  dock  drift.  The  AFGL  stations  were  also 
connected  by  radio  antenna  to  WWVB. 

Deployment  of  field  stations  consisted  of  burying  the  seismometer  approximately  one 
foot  below  the  ground  surface.  The  sensor  was  oriented  so  that  the  "Longitudinar  arrow  pointed 
towards  Magnetic  North.  The  sensor  was  then  leveled  and  connected  to  one  of  the  digital 
recorders.  A  WWVB  radio  antenna  was  also  connected  to  the  recorder  as  a  back-up  time  source. 
Care  was  taken  to  set  up  the  same  sensor  and  recorder  pair  at  the  same  site  and  calibrations 
were  made  at  each  site.  At  the  scheduled  blast  time,  field  recorders  were  switched  on  manually 
to  record  the  event.  Stations  were  manually  run  because  most  of  them  would  either  not 
automatically  trigger  on  an  event  or  else  the  site  would  be  so  "noisy"  that  too  many  folse 
triggers  would  use  up  all  of  the  tape.  Also,  blasting  times  tend  to  be  too  unpredictable  to  leave 
the  recorders  on  trigger  mode.  Therefore  after  the  recorders  were  turned  on,  the  person 
monitoring  the  site  would  call  the  quarry  scale  house  to  confirm  whether  or  not  the  blast  hsKl 
taken  place. 
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1 .5  Description  of  Data  Set 

The  Actors  that  affect  the  quality  and  success  with  which  a  seismic  event  is  recorded  are 
innumerabie.  To  discuss  every  one  of  those  factors  in  detail  would  stray  from  the  focus  of  this 
thesis.  However,  on  a  most  fundamental  level,  an  investigation  of  this  type  requires  that  at  least 
the  following  three  criteria  be  satisfied  for  each  seismogram:  1)  accurate  timing;  2)  a 
relatively  high  signal-to*noise  ratio;  and  3)  precise  source  and  receiver  locations.  As 
straightforward  as  that  may  sound,  obtaining  all  three  of  those  criteria  for  a  single  seismogram 
can  be  somewhat  of  a  challenge. 

The  original  data  set  used  to  begin  work  consisted  of  approximately  110  seismograms 
that  crossed  about  49  different  paths  throughout  the  study  area.  Several  paths  that  were  part  of 
the  1989  summer  field  experiment  had  to  be  eliminated  immediately  because  of  internal  dock 
failures.  For  this  type  of  investigation,  data  without  accurate  timing  is  useless.  The  remaining 
seismograms  were  plotted  using  a  band-pass  filter  with  a  high-pass  cutoff  of  12.5  Hz  and  a 
low-pass  cutoff  of  .01  Hz  to  filter  out  some  of  the  badcground  noise  and  also  to  enhance  the  event 
signal.  At  this  point  four  more  paths  were  eliminated  because  these  (MIT)  network  stations  did 
not  adequately  record  the  event.  The  reason  for  this  could  have  been  that  the  distance  between 
source  and  receiver  was  too  great,  and/or  that  the  signal-to-noise  ratio  was  so  low  that  the 
event  could  not  be  distinguished  from  background  noise. 

After  analyzing  the  waveforms  on  the  remaining  seismograms,  the  data  set  was  divided 
into  the  following  four  categories  according  to  the  quality  of  the  Rg  signals; 

Category  1)  Excellent  quality.  Seismogram  has  excellent  Rg  and  P  and/or  S 
waveforms.  High  signal-to-noise  ratio. 

Category  2)  Good  quality.  Seismogram  has  good  Rg  waveform,  but  may  have  a 

significant  amount  of  background  noise  that  interferes  with  the  signal. 

Category  3)  Fair  quality.  Seismogram  has  average  to  poor  waveform  and  a  low 
signal-to-noise  ratio. 
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Category  4)  Poor  quality.  Seismogram  has  no  apparent  Rg  wave.  Other  waveforms 
are  of  poor  quality  or  nonexistent.  Difficult  to  distinguish  signal  from 


background  noise. 

Table  2  is  a  list  of  all  the  data  files  and  their  corresponding  categories.  When  processing 
for  group  velocity,  it  was  found  that  those  seismograms  in  category  4  did  not  yield  useful 
measurements.  Therefore,  those  paths  were  eliminated  from  the  Rg  dispersion  analysis.  What 
remained  after  all  these  inspections  was  a  data  set  consisting  of  approximately  81  seismograms 
that  trace  44  different  paths  within  the  study  area.  Of  the  81  seismograms,  17  are  quality  1; 
22  are  quality  2;  28  are  quaiity  3;  and  14  are  quality  4  seismograms  Figures  7  -  10  are 
examples  of  seismograms  that  are  representative  of  their  respective  categories.  The  paths 
range  in  length  from  about  2  km  to  about  82  km.  For  a  few  of  the  paths  there  are  more  than  one 
recording.  In  such  cases  an  analysis  was  done  on  the  multiple  recordings  separately.  Multiple 
data  was  then  averaged  for  a  given  path  and  examined  that  way  as  well.  In  the  final  analysis,  it 
was  elected  to  display  the  data  as  separate  data  files. 
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Table  2 


Station 

Category 

Station 

Category 

LC1 

2 

87202-GLO 

1 

LC2 

1 

87210-ONH 

3 

LC3 

1 

87210-WFM 

3 

LjC4 

1 

87210-WNH 

4 

AZ1 

2 

87210-DNH 

3 

AZ2 

1 

87210-GLO 

1 

AZ3 

1 

88089-ONH 

4 

AZ4 

1 

88089-PNH 

3 

NSAO 

2 

88089-WFM 

2 

NSA1 

2 

88089-ONH 

3 

NSA2 

2 

88089-GLO 

2 

NSA3 

2 

88089-UXB 

2 

NSB1 

2 

89177-ONH 

4 

NSB2 

2 

89177-PNH 

4 

NSB3 

2 

89177-WFM 

3 

NSB4 

2 

89177-DNH 

2 

NSB5 

2 

89177-GLO 

1 

Plot 

4 

89192-ONH 

4 

P102 

3 

89192-PNH 

2 

PI  03 

3 

89192-WFM 

2 

P104 

4 

89192-DNH 

3 

P107 

3 

89192-GLO 

3 

P1 10 

1 

89192-NW01 

1 

P112 

1 

89192-NW02 

3 

P113 

2 

89192-NW06 

2 

P115 

1 

89192-NW10 

1 

CLA 

1 

89208-NW03 

3 

WES 

1 

89208-NW04 

3 

P201 

4 

89208-NW05 

3 

P202 

4 

89213-DNH 

4 

P203 

4 

89213-GLO 

3 

P204 

4 

89213-WFM 

2 

P205 

3 

89213-NW01 

3 

P206 

3 

89213-NW03 

3 

P301 

4 

89213-NW04 

2 

P304 

3 

89213-NW05 

3 

P305 

3 

89213-NW07 

2 

87202-ONH 

4 

89213-NW08 

3 

87202-PNH 

3 

89213-NW10 

3 

87202-WFM 

87202-DNH 

1 

3 

89213-NW11 

3 
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Z  GBXOGYAND  CRUSTALSTRUCTURE  0FSCXJ7VCRN 
2.1  Geology  of  Southern  New  England 

Skehan  and  Osberg  (1979)  descdbed  the  Northern  Appalachian  Mountain  Region  as 
geologically  complex  and  comprised  of  regions  with  distinct  stratigraphy,  structure  and 
deformation.  They  suggested  that  this  region  is  a  series  of  *lithotectonic  terranes"  that  have 
become  sutured  together  as  a  result  of  geologic  activity  during  the  Middle  Paleozoic  through 
Mesozoic  times.  The  terranes,  which  are  believed  to  at  one  time  have  been  separate  and 
individual  units,  came  together  as  a  result  of  massive  plate  collisions.  This  created  a  "collage" 
of  geologic  subdivisions,  each  with  distinct  characteristics. 

Figure  11  shows  the  terranes  of  SNE.  They  are  from  west  to  east:  (1)  the  Proto-North 
American  terrane  (NA)  which  represents  Late  Proterozoic  and  Early  Paleozoic  rift  and  shelf 
sedimentary  rocks  deposited  on  Grenville  basement;  (2)  the  Connecticut  Valley  (CV)  sequence  of 
thrust  slices  separated  from  NA  by  Cameron's  line  (CL);  (3)  the  Hartford  Rift  Basin  (HRB) 
which  consists  of  Mesozoic  red  bed  sediments  and  basaltic  igneous  rocks  covered  by  relatively 
thick  Pleistocene  glacial  sediments;  (4)  the  Bronson  Hill  (BH)  Ordovician  volcanic  arc;  (5)  the 
Kearsarge-Central  Maine  Synclinorium  (KC),  dominantly  Silurian-Devonian  metasedimentary 
rocks  intruded  by  Devonian  piutons;  (6)  the  Merrimac  Trough  (MT),  metasedimentary  rocks  of 
uncertain  age  separated  from  adjacent  terranes  by  faults  ;  (7)  the  Putnam-Nashoba  terrane 
(PN)  also  bounded  by  faults  and  consisting  in  part  of  Late  Proterozoic  age  rocks;  and  (8)  the 
Avaionian  terrane  (AV),  Late  Proterozoic  volcanic  rocks  with  a  cover  of  distinctive 
fossiliferous  sedimentary  rocks  (Skehan  and  Kafka,  1990). 

The  Bronson  Avalon  Dispersion  Region  (BADR)  encompasses  several  of  these 
subdivisions  in  that  it  extends  from  the  Bronson  Hill  Anticlinorium  (BH)  eastward  to  the 
Avalonian  Terrane  (AV).  Within  this  area,  the  basement  rocks  are  characterized  by  Late 
Precambrian  gneisses  and  quartzose  sedimentary  rocks  with  layered  gneisses  and  plutonic 
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NA  =  Proto-North  American  Terrene  KC  =  Kearsarge-Central  ME 
CV  =  Connecticut  Valley  Synclinorium  Syncllnoiium 

CL  =  Cameron's  Line  MT  =  Merrimack  Trough 

HRB  =  Hartford  Rift  Basin  PN  =  Putnam-Nashoba  Terrene 

BH  =  Bronson  Hill  Anticlinorium  AV  =  Avalonian  Superterrane 


Figure  11:  Tectonic  terrenes  of  Southern  New  England,  (from  Kafka 
and  Skehan,  1990) 
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gneisses  that  are  partly  Precambrian  and  partly  Ordovician  in  age.  Late  Precambrian  rocks  are 
exposed  in  the  cores  of  the  Pelham  dome,  Clinton,  Stony  Creek,  and  Lyme  domes  and  also  in  the 
Willimantic  Dome  of  southeastern  CT  arxt  the  Milford  Antidinorium  of  eastern  MA  and  Rl  (Hall 
and  Robinson,  1982). 

Middle  Ordovician  strata  within  tfte  BADR  lie  directly  on  dome  gneisses  and  consist 
mainly  of  suKidic  mica  schist,  although  there  are  rocks  of  volcanic  derivation  including 
amphibolites  and  felsic  gneisses  that  s^e  interspersed  with  mica  schists.  The  Middle  Ordovician 
schists  are  mostly  exposed  east  of  the  Bronson  Hill  Antidinorium.  Elsewhere  the  age  of  similar 
rocks  is  less  certain.  These  rocks  indude  metamorphosed  mafic  and  felsic  voicanics  such  as  the 
Marlboro  Formation  in  MA  and  the  Quinebaug  in  CT  that  have  been  dated  by  Hill  and  others 
(1984)  as  in  the  range  550  to  450  million  years  old.  The  voicanics  are  in  contrast  with  and 
possibly  overlain  by  thick  units  of  peltic  schists  and  gneisses  such  as  the  Nashoba  Formation  in 
MA  and  the  Tatnic  Hill  Formation  in  CT  (Hall  and  Robinson,  1982). 

Silurian  rocks  of  the  BADR  consist  of  a  quartz  pebble  type  of  conglomerate-quartzite 
called  the  Clough  Quartzite,  in  northern  MA  this  rock  contains  Llandovery  fossils.  The  Fitch 
Formation  that  extends  into  northwestern  NH  also  contains  fossils.  These  fossils  and  the  tens- 
like  distrfoution  of  the  rock  suggest  a  marine  or  shoreline  type  depositional  environment. 
Moving  eastward,  the  Silurian  rocks  have  been  interpreted  as  thickening  into  a  section  of 
feldspathic  calcareous  granulite  sulfkfic  schist  (i.e.  the  Paxton  Formation).  The  eastern  limit  of 
this  Silurian  basin  exposure  is  along  foe  Clinton-Newbury  Fault  in  MA  and  the  top  of  the  Tatnic 
Hill  Formation  in  CT  (Hall  and  Robinson,  1982). 

Devonian  rocks  in  foe  BADR  are  mostly  schist  interlayered  with  quartzose  marble.  The 
western  BADR  is  dominated  by  gabbro  to  granite  plutons  such  as  the  Coys-Hill  Cardigan  Pluton, 
the  Ashuelot  Pluton,  and  parts  of  foe  Fitchburg  Pluton.  The  eastern  BADR  rocks  are 
predominantly  gray-weathered  graphite  schists  of  the  Littleton  Formation  (Hall  and  Robinson, 
1982). 
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Carboniferous  secfimentary  rocks  are  exposed  in  three  localities  west  of  the  Clink>n> 
Newbury  Fault  in  Central  MA  (Skehan,  Rast.  and  Mosher  1986).  Carboniferous  to  Permian 
plutons  include  the  Pinewood  Adamellite  Piuton  in  southeastern  CT,  the  Narragansett  Pier 
Granite  of  Rl.  and  the  Milford  Granite  in  southern  NH  (Hall  and  Robinson.  1982). 

Within  the  BADR  there  are  four  structural  subdvisions.  They  are  from  west  to  east:  the 
Bronson  Hill  Antidinorium.  the  Kearsarge-Central  Maine  Synciinorium,  the  Merrimack 
Trough,  the  Nashoba  Terrane.  and  the  Avalon  Terrane.  Within  the  synciinorium.  the  dominant 
structures  are  upright  isoclinal  folds,  some  of  which  have  been  deformed.  Evidence  of  inverted 
rock  sequences  suggest  that  there  has  been  some  recumbent  folding.  The  Bronson-Hili 
Antidinorium  is  a  complex  structure  distorted  by  domes  and  basins  and  is  cut  along  the  west 
margin  by  Late  Paleozoic  and  Mesozoic  faults  (Skehan  and  Osberg.  1979). 

The  BADR  is  also  distinguished  by  many  well  documented  faults.  The  Ciinton-Newbuty 
Fault  traces  southwesterly  Into  the  Lake  Char  Fault  near  Worcester.  MA.  It  can  then  be  traced 
into  the  Honey  Hill  Fault  of  CT.  Four  faults  have  been  mapped  that  splay  off  to  the  north  of  the 
Clinton  Newbury  Fault  into  NH.  They  are.  from  west  to  east,  the  Pinnacle.  Hall  Mountain* 
Campbell  Hill.  Silver  Lake  and  RInt  HHI  Faults.  The  Flint  Hill  and  Silver  Lake  Faults  merge  to 
the  south  to  become  the  Wekepeke  Fault  in  MA.  Eventually  the  Wekepeke  Fault  joins  the  Clinton- 
Newbury  Fault  near  Worcester.  MA  (Lyons  et  ai.  1982). 

2.2  Crustal  Structure  of  Southern  New  England 

Studies  of  the  crustal  velocity  structure  in  the  New  England  area  are  numerous.  Past 
work  using  body  wave  or  surface  wave  methods  to  determine  seismic  velocity  structure  indude 
for  example  Linehan  (1962).  Breitling  (1965).  Chiburis  et  al.  (1977).  Bruenig  (1980). 
and  Taylor  and  Toksoz  (1982).  More  recently,  studies  using  group  velodty  dispersion  of  Rg 
waves  to  constrain  variations  in  the  seismic  structure  of  the  shallow  crust  include  Kafka  and 
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Dollin  (1985).  Kafka  (1985),  Kafka  and  McTigue  (1985),  McTigue  (1986),  Saikia  et  al. 
(1986),  Gnewuch  (1987),  Kafka  (1988),  and  Tu  (1989). 

Kafka  and  Dollin  (1985)  studied  group  velocity  dispersion  for  Rg  waves  with  periods 
between  0.5  and  2.0  sec  that  were  recorded  from  quarry  blasts  in  SNE.  In  that  study,  they 
identified  five  cSstinct  dispersion  regions.  Normal  dispersion  was  observed  in  the  period  range 
of  0.5-1 .5  sec  indicating  that  there  must  be  a  reiatively  low  velocity  surfidal  iayer  (about  0.5 
to  1 .0  km  thick)  in  the  upper  crust  throughout  the  study  area.  The  seismic  velocities  in  this 
ayer  displayed  some  lateral  variation  and  appeared  to  oorreiate,  to  some  extent,  with  the 
geoiogy. 

McTigue  (1986)  used  groi^  velodty  dispersion  of  Rg  waves  between  periods  of  0.5  and 
2.0  sec  to  map  lateral  and  verticai  variations  in  the  seismic  velocity  structure  of  the  crust 
beneath  SNE.  He  noted  the  same  five  (fispersion  regions  as  Kafka  and  Skehan  (1990)  and 
obtained  velocity  models  approximating  the  upper  crustal  structure  for  each  region.  All  Rg 
waves  analyzed  in  his  study  displayed  normal  dispersion  at  shorter  periods,  again  indicating  the 
presence  of  a  low  velocity  layer  about  0.5  to  1.0  km  thick  near  the  surface.  He  concluded  that 
the  Avalonian  Terrane  had  a  slightly  different  crustal  structure  than  the  area  to  the  west.  His 
velocity  models  suggested  to  him  that  the  Avalonian  extends  beneath  the  Nashoba  Terrane  up  to 
50  km  west  of  the  westernmost  exposure  of  the  Avalonian  terrane.  McTigue  concluded  that  this 
may  represent  a  subsurface  extension  of  the  Douglas  Woods  Antidlne. 

Saikia  et  al.  (1986)  measured  groi^  velocities  from  New  England  quarry  blasts  along 
various  paths  with  distances  up  to  180  km  from  the  source.  They  inverted  observed  dispersion 
data  to  obtain  shear  wave  velocity  models  for  various  areas  in  SNE.  They  also  computed 
theoretical  dispersion  curves  and  displacement  depth  eigenfunctions  for  Rg  waves  using  various 
earth  models.  Synthetic  Rg  waves  were  ttien  calculated  for  explosion  sources  using  those  earth 
models  from  which  the  theoretical  dispersion  cunres  and  the  observed  dispersion  curves  had  the 
best  fit. 
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Farther  north,  in  southeastern  Maine,  Kafka  and  Reiter  (1987)  analyzed  the  dispersive 
properties  of  Rg  waves  that  were  recorded  from  blasts  detonated  by  the  U.S.  Geological  Survey  as 
part  cf  the  Maine  Seismic  Refraction  Experiment  in  1984.  Group  velocities  for  the  period 
range  of  0.4  to  1 .6  sec  are  sensitive  to  the  velocity  structure  at  depths  ranging  from  near  the 
surface  to  a  depth  of  approximately  2  km.  The  group  velocities  that  were  observed  in  this  study 
were  found  to  be  different  from  one  path  to  the  next  The  foct  that  the  group  velocities  were 
significantly  higher  for  paths  parallel  to  the  grain  of  the  Appalachians  than  for  paths  transverse 
to  the  grain,  suggested  to  them  that  the  shallow  crustal  structure  beneath  the  Appalachians  of 
southeastern  Maine  is  laterally  anisotropic.  Rg  group  velocities  were  inverted  to  obtain  shear 
wave  velocity  models  of  the  upper  2.0  km  of  the  crust.  A  Vp  model  was  constructed  using  a  Vp/Vs 
ratio  of  1.78.  Their  resulting  models  indicated  that,  on  average,  Vp  ranged  from  about  4.9 
km/sec  in  the  upper  0.3  km  to  6.2  km/sec  at  a  depth  of  about  2  km.  The  obsen/ed  range  of  Rg 
groi4)  velocities  in  this  study  is  evidence  that  there  is  a  variation  in  shallow  crustal  structure 
that  is  dependent  on  either  path  location,  path  orientation,  or  both. 

Kafka  (1988)  summarized  Rg  dispersion  studies  in  SNE  and  concluded  that  although  Rg 
dispersion  varies  significantly  throughout  SNE,  there  is  a  rather  large  region  where  the  Rg 
dispersion  is  curiously  uniform.  That  region  is  the  BADR,  which  was  given  that  name  based  on 
the  findings  in  those  studies.  The  Rg  observations  suggest  that  on  the  scale  of  variations 
deDneated  by  these  Rg  studies,  the  shallow  crust  beneath  this  area  is  homogeneous.  The 
observed  variations  in  the  shallow  crustal  structure  across  SNE  is  summarized  by  Kafka  and 
Skehan  (1990). 

Tu  (1989)  studied  quarry  blasts  in  New  Hampshire  and  Vermont  that  were  recorded  by 
the  NESN,  and  determined  group  velocities  for  Rg  waves  with  periods  between  0.3  sec  and  2.0 
sec.  The  paths  used  in  this  study  ranged  from  about  10  to  168  km  in  length.  Results  of  this 
study  show  that  group  velocities  in  this  area  of  New  Hampshire  and  Vermont  very  much 
resemble  those  of  the  BADR,  which  suggests  that  the  BADR  extends  into  New  Hampshire  and 
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adjacent  Vermont.  However,  because  Kaika  and  Reiter  (1987)  discovered  significant  variation 
in  seismic  veiodty  structure  in  southeastern  Maine,  there  must  be  a  "transition  zone* 
somewhere  between  southern  New  Hampshire  and  southeastern  Maine  where  significant  iateral 
variation  begins  to  occur. 
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3.  SURFACE  WAVE  THEORY 
3.1  Rg  Waves 

An  Rg  wave  is  a  short-period  surface  wave.  Surtace  waves  are  a  type  of  seismic  wave 
that  traveis  along  the  free  surface  of  the  earth.  A  Rayleigh  wave  is  a  type  of  a  surtace  vrave  in 
which  the  ground  displacement  is  in  the  vertical  plane,  with  a  retrogr^Je  ell^tical  partide 
motion.  Fundamental  mode  Rayleigh  %vaves  with  periods  between  about  0.2  sec  and  2.2  sec  (Rg) 
are  very  often  observed  on  seismograms  of  quarry  blasts  and  shallow-focus  earthc^akes. 
Quarry  blasts  and  other  explosions  that  are  detonated  at  or  just  below  the  earth's  surface 
usually  generate  strong  Rg  wave  signals,  making  them  ideal  seismic  sources  for  the  study  of  Rg 
wave  dispersion.  The  dispersive  properties  of  an  Rg  wave  are  sensitive  to  lateral  variations  in 
the  upper  few  kilometers  of  the  crust. 

Because  the  major  component  of  Rayleigh  wave  partide  motion  is  shear  in  the  vertical 
plane  and  the  minor  component  is  oompresdonal-dilitational  in  the  horizontal  plane,  Rayleigh 
waves  are  most  sensitive  to  variations  in  shear  wave  velodty  (Dobrin  et  al,  1951).  Panza  and 
Calagnile  (1975)  calculated  theoretical  phase  and  group  velocity  dispersion  curves  for 
continental  earth  structures,  and  they  demonstrated  that  Rg  can  be  modeled  as  a  fundamental 
mode  Rayleigh  wave.  Displacement-depth  eigenfunctions  of  fundamental  mode  Rayleigh  waves 
computed  for  the  Chburis  (1977)  crustal  model  in  the  period  range  of  0.2  -  2.5  sec  show  that 
Rg  wave  particle  displacement  is  at  a  maximum  near  the  surface  and  decreases  to  essentiatiy 
zero  at  a  depth  of  about  10  km.  it,  therefore,  makes  sense  that  surface  explosions  such  as 
quarry  blasts  and  shallow-focus  earthquakes  are  excellent  sources  for  the  generation  of  surfoce 
waves  (see  Kafka,  1990). 

Observations  of  Rg  wave  propagation  have  shown  that  different  frequency  components  of 
the  wave  travel  with  different  velocities.  This  phenomenon  causes  the  shape  of  the  wavetrain  to 
change  with  distance  from  the  source.  When  this  happens,  the  wave  is  said  to  be  dispersive 


28 


(Dobrin,  1951).  Because  Rg  waves  are  dispersive,  two  kinds  of  measurements  for  veiocity  can 
be  made.  The  phase  vetodty  of  an  Rg  wave  is  the  measurement  of  the  distance  traveled  per  unit 
of  time  by  a  point  of  constant  phase,  such  as  a  peak  or  a  trough.  If  you  were  to  draw  an  envelope 
around  the  peak  of  a  wave  pulse  and  measure  the  speed  with  which  the  wave  packet  was 
traveling,  then  you  would  be  measuring  the  group  velocity  of  that  wave.  The  group  veiocity  of  a 
wave  can  be  defined  as: 

U  =  C-X^  =  C+(o|5- 

dcD 

where  C  =  phase  velocity.  X  =  wavelength,  and  cd  s  angular  frequency. 

3.2  Measuring  Group  Velocities 

Group  velocities  of  Rg  waves  are  measured  using  a  technique  of  narrow  band  pass 
filtering  developed  by  Dziewonski  et  ai  (1969).  This  method,  which  has  been  frequently  used 
and  upgraded  by  seismologists  at  Weston  Observatory,  is  used  to  study  the  variation  of  wave 
energy  as  a  function  of  velocity  and  period.  The  signal  is  first  Fourier  transformed  and  then  a 
Gaussian  filter  is  applied  in  narrow  frequency  bands  at  a  series  of  center  frequencies.  For  each 
center  frequency,  the  group  arrival  time  is  approximated  by  measuring  the  arrival  time  of  the 
maximum  amplitude  of  the  envelope  of  the  filtered  signal.  This  process  is  repeated  for  a  series 
of  center  frequencies  and  group  velocities  are  estimated  at  each  center  frequency. 

The  result  of  this  analysis  is  a  contour  plot  of  group  velocity  vs.  period  showing  the  log 
(base  10)  of  energy  for  each  center  frequency  as  a  function  of  group  vefocrty  (Figure  12).  The 
shape  of  the  contour  plot  corresponds  to  the  shape  of  the  energy  envelope  at  each  center 
frequency  (Dollin,  1984).  The  group  velocity  curve  that  results  from  the  plot  represents  the 
velocities  at  which  maximum  energy  arrived  for  each  center  frequency. 

The  example  shown  in  Rgure  12  represents  the  contour  plot  for  a  blast  that  was 
detonated  at  the  San-Vel  quarry  and  recorded  at  Weston  Observatory's  network  station  (WES). 
The  lower  numbers  on  the  plot  represent  higher  amplitudes. 
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Figure  ^2  ;  Contour  plot  of  San-Vel  quarry  blast  recorded  at 
station  WES  with  group  velocity  curve  superimposed-  Below  is 
seismogram  of  the  event. 
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4.  DATAPROCeSSINQ 
4.1  Origin  Times 

in  order  to  analyze  group  velocities  of  Rg  waves,  accurate  estimates  of  event  origin  times 
are  critical.  Fortunately,  for  this  investigation  ail  of  the  data  gathered  prior  to  1989  had  very 
accurate  origin  time  estimates.  For  those  events,  an  accelerometer  was  placed  in  the  quarry  as 
dose  to  the  shot  point  as  possible  (at  a  dtetance  of  approximately  100  yards).  The 
accelerometer,  which  was  connected  to  a  digital  recorder  like  all  of  the  other  field  sites, 
recorded  the  shot  on  tape.  Because  of  the  dose  proximity  of  the  sensor  to  the  event  we  were 
able  to  estimate  the  origin  time  to  within  a  few  hundredths  of  a  second  by  measuring  the  arrival 
time  of  the  first  motion  on  the  record.  Table  3  lists  the  1987>1988  events  with  their  origin 
times. 

Origin  times  for  the  1989  field  data  and  events  thereafter  were  estimated  to  within 
about  0.1  to  0.2  sec  of  the  actual  origin  time  using  a  statistical  method  based  on  all  of  the  data 
from  events  with  known  origin  times.  To  do  this,  travel  times  (TT)  were  calculated  for  all  pre- 
1989  seismograms  by  subtracting  the  event  origin  time  (OT)  from  the  P-wave  arrivai  time 
(P-AT).  These  calculations  were  entered  into  the  StatView  SE-i-Qraphics  statistical  software 
package  designed  by  Abacus  Concepts,  Inc.  for  the  Madntosh  computer.  Using  this  program,  the 
travel  times  were  plotted  as  a  function  of  distance  (Figure  13).  A  linear  regression  performed 
on  these  data  yieided  an  equation  for  travel  time: 

TT-.171(x)  +  .151,  (2) 

where  x  >  distance  from  source  to  receiver.  Using  this  equation,  theoretical  travel  times  could 
be  calculated  for  ail  of  the  seismograms  with  unknown  origin  times.  Origin  times  for  these 
events  were  then  estimated  by  subtracting  the  theoretical  travel  time  from  the  actual  P-wave 
arrival  time  (picked  directly  off  of  the  seismogram).  The  final  estimated  origin  time  used  for 
the  processing  of  this  data  was  determined  by  averaging  all  of  individual  origin  time  estimates 
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Table  3 

EVENTS  WITH  KNOWN  ORIGIN  TIME 

EVENTLOCATION  LAT  LON 

SAN-VEL/LONESTAR  42.554N  -71.517W 

QUARRY 


YEAR  JUUAN  DAY  ORIGIN  TIME 


1987 

202 

16:58:14.97 

1987 

210 

18:13:56.17 

1988 

089 

18:01:50.31 
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TABLE 4 


ESTIMATING  ORIGIN  TIMES 

USING  THE  EQUATION:  TT  =  .171  (X)  +  .151  (from  figure  4.1.1) 

JULIAN  DAY: 

STATION  DIST.  (km^  P-AT _ H _ OUSSL^ 


DAY  177 


NW02 

7.86 

18:28:56.72 

1.50 

18:28:55.22 

WFM 

6.65 

18:28:56.54 

1.29 

18:28:55.25 

WES 

24.71 

18:28:59.84 

4.38 

18:28:55.46 

Average  estimated  origin  time: 

18:28:55.36 

DAY  192 
NWOl 

5.87 

19:33:55.88 

1.15 

19:33:54:73 

NW02 

7.86 

19:33:56.30 

1.50 

19:33:54.80 

NW06 

17.64 

19:33:58:05 

3.17 

19:33:54.88 

Average  estimated  origin  time: 

19:33:54.80 

DAYJOa 

NW03 

8.96 

17:29:49.83 

1.68 

17:29:48.15 

NW04 

13.80 

17:29:50.60 

2.51 

17:29:48:09 

NW05 

15.69 

17:29:51.12 

2.83 

17:29:48.29 

Average  estimated  origin  time: 

17:29:48:18 

DAY  213 
NWOl 

5.87 

19:29:26.91 

1.15 

19:29:25.76 

NW03 

8.96 

19:29:27.52 

1.68 

19:29:25.84 

NW04 

13.80 

19:29:28.26 

2.51 

19:29:25.75 

NW05 

15.69 

19:29:28.88 

2.83 

19:29:26.05 

NW07 

20.25 

19:29:29.62 

3.61 

19:29:26.01 

NW08 

22.79 

19:29:29.98 

4.05 

19:29:25.93 

Average  estimated  origm  time:  19:29:25.89 
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Table  5 


EVENT  (YRDAY)  ORIGIN  TIME 


87202 

16:58:14.97 

87210 

18:13:56.17 

88089 

18:01:50.31 

89177 

18:28:55.31 

89192 

19:33:54.80 

89208 

17:29:48.18 

89213 

19:29:25.89 

Table  6 


Event 

U 

P-vel  ttheo 

MI 

(sec) 

(km/sec) 

(sec) 

(sec) 

89177-NW02 

1.36 

5.78 

1.50 

.14 

89192-NW01 

1.08 

5.44 

1.15 

.07 

89192-NW02 

1.51 

5.21 

1.50 

.01 

89192-NW06 

3.25 

5.43 

3.17 

.08 

89208-NW03 

1.65 

5.43 

1.68 

.03 

89208-NW04 

2.42 

5.70 

2.51 

.09 

89208-NW05 

2.94 

5.34 

2.83 

.11 

89213-NW01 

1.02 

5.76 

1.15 

.13 

89213-NW03 

1.63 

5.50 

1.68 

.05 

89213-NW04 

2.37 

5.82 

2.51 

.14 

89213-NW07 

3.73 

5.43 

3.61 

.12 

89213-NW08 

4.09 

5.57 

4.05 

.04 
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DISTANCE 


Figure  13  :  Travel  time  plotted  as  a  function  of  distance. 
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for  each  event.  Table  4  shows  examples  of  events  where  origin  time  was  estimated  using  this 
procedure. 

Using  this  method,  origin  times  were  estimated  for  aH  of  the  suit  ner  field  events.  Table 
5  is  a  list  of  all  seven  events  recorded  in  the  field  between  1987  and  1989  with  their 
respective  origin  times.  The  effect  of  origin  time  errors  on  group  veiocity  measurements  is 
discussed  in  Section  42. 

When  using  a  statistical  method  to  estimate  any  measurement  there  is  bound  to  be  some 
margin  of  error.  As  mentioned  above,  accurate  origin  times  are  criticai  when  processing  a 
seismogram  for  group  velocity.  One  thing  that  must  be  considered  is,  ”What  is  the  margin  of 
error  in  estimating  origin  times?"  To  test  the  accuracy  of  the  statistical  method  described 
above,  the  same  method  was  applied  to  estimate  the  origin  times  for  events  with  known  origin 
times.  When  these  estimated  origin  times  were  compared  with  the  actual  origin  times,  the 
method  was  found  to  be  accurate  to  within  .1  to  2  sec.  whereas  the  1987  and  1988  events  are 
probably  accurate  to  within  a  few  hundredths  of  a  second.  Also.  H  you  calculate  travel  times  for 
the  data  using  the  estimated  origin  times  and  the  actual  origin  times  and  compare  the  travel 
times  for  both  groups,  you  can  see  that  the  difference  in  travel  time  ( A  tt)  for  each  event  is 
within  +/-  .14  sec.  Table  6  lists  the  events,  their  actual  travel  times  (tt),  P-wave  velocity 
(P-vel),  their  theoretical  travel  times  (ttheo)  and  the  difference  between  the  two  (Att). 

42  Effect  of  Origin  Time  Errors  On  Group  Velocity  Measurements 

It  is  important  to  note  that  errors  in  origin  time  are  going  to  have  a  greater  effect  on 
paths  of  shorter  length  than  on  those  that  are  longer.  This  is  demonstrated  in  the  following 
example.  Consider  a  station  whose  distance  from  the  source  (x)  is  30  km.  Assuming  an  Rg  wave 
velocity  (U)  of  2.5  km/sec,  the  Rg  wave  travel  time  would  be  12  sec.  if  we  add  an  error  of  0.1 
sec  to  the  arrival  time,  we  would  measure  a  travel  time  of  12.1  sec.  This  in  turn  would  yield  an 
Rg  veiocity  of  U'  *  2.48  km/sec.  The  differences  between  the  two  velocity  measurements  would 
be  A  U  =  0.02  km/sec.  If  we  set  x  »  60  km.  then  U*  «  2.49  km/sec  and  A  U  would  be  equal  to 
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0.01  km/sec.  For  the  same  0.1  sec  error,  the  differences  or  errors  In  velocity  measurements 
are  cut  in  half  as  we  increased  the  distance  from  source  to  receiver  by  two.  If  we  continued  to 
increase  the  distance  (x),  then  the  errors  in  velocity  would  become  less  and  less  significant. 
Conversely,  as  distance  decreases  the  origin  time  errors  have  a  greater  effect  on  the  accuracy  of 
the  group  velocity  measurements.  The  graph  in  Figure  14  shows  the  change  in  group  velocity 
with  increasing  distance  for  three  different  error  measurements. 

As  discussed  above,  timing  is  essential  to  the  study  of  seismic  velocities.  Great 
care  is  taken  to  synchronize  very  sensitive  time  keeping  devices  in  all  instruments  used 
in  a  seismic  experiment.  Curiously,  when  processing  this  data  set  some  very  suspicious 
group  velocity  measurements  were  found  from  what  appear  to  be  some  very  good 
seismograms.  Because  there  was  no  obvious  explanation  for  these  suspicious 
measurements,  the  first  step  I  taken  to  investigate  this  problem  was  to  ensure  that  the 
file  was  processed  correctly.  In  doing  so,  a  check  was  made  of  the  timing  of  the  event. 

The  NESN  and  MIT  networks  are  set  up  so  that  every  time  an  event  is  recorded,  a  time 
file  is  generated  that  contains  the  time  foat  is  read  off  of  a  satellite  dock  and  is  then 
correlated  fo  the  event  trace.  The  program  that  generates  the  time  file  calculates  at  what 
time  the  first  data  point  arrives  and  lines  the  trace  up  with  the  time  from  there,  if  the 
program  does  not  accurately  calculate  the  time  of  the  first  data  point  on  the  trace,  then 
the  timing  for  the  entire  event  will  be  off.  When  I  plotted  the  time  files  for  the  above 
mentioned  events.  I  found  that  indeed  the  second  and  minute  marks  did  not  line  up 
precisely,  as  they  should.  Each  event  was  offset  by  a  small  amount  which  I  have  noted  in 
Table  7.  I  processed  each  file  again,  having  corrected  for  the  time  discrepandes.  As  a 
result  of  the  corrections,  each  seismogram  has  yielded  usable  data. 
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Rgure  14:  Change  in  group  velocity  with  increasing 
distance  for  examples  of  different  origin  time  errors. 


Table? 


^ar _ Qay _ Clock  Error _ Station  Affiliation 


1987 

202 

.04  sec  (slow) 

MIT 

1987 

210 

.15  sec  (slow) 

MIT 

1988 

089 

.71  sec  (fast) 

MIT 

1989 

177 

.30  sec  (slow) 

MIT 

1989 

192 

.30  sec  (slow) 

MIT 

1989 

213 

.30  sec  (slow) 

MIT 
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5.  RESULTS 

5.1  Summary  of  Results  from  Previous  Studtos  of  SNE 

As  a  result  of  ttiis  study,  observed  Rg  group  vetoc^ies  in  eastern  MA  and  southern  NH 
range  between  2.00  to  3.29  km/sec  for  a  period  range  of  0.2  to  2.2  sec.  Compared  with  the 
results  of  previous  investigations  of  Rg  wave  group  velocity  dtepersion  (T able  8).  it  was  found 
that  the  range  of  group  velocities  from  this  study  is  most  simiiar  to  that  observed  by  Kaflta  and 
Dollin  (1985)  for  paths  in  Southern  New  England  (2.0  to  3.3  km/8ec). 

Originally,  the  results  of  Kafka  and  DoIRn  (1985)  revealed  lateral  variations  in  group 
velocity  dispersion  cun/es  that  scattered  around  three  distinguishable  ranges  of  group 
velocities.  Group  velocities  for  paths  that  crossed  southwestern  Connecticut  were  noticeably 
higher  than  those  paths  that  traversed  eastern  Connecticut.  Both  of  those  groiq)s  of  paths  were 
foster  than  paths  that  crossed  the  Hartford  Rift  Basin.  As  a  result,  Kafka  &  Dollin  designated 
each  area  as  a  distinct  dispersion  region.  More  recently,  however,  Kafka  and  Bowers  (1991) 
discovered  that  in  the  Kafka  &  Dollin  (1985)  study  an  incorrect  source  to  receiver  distance  was 
applied  in  processing  paths  that  extended  from  the  various  quarries  to  statfon  BCT  (Brookfield, 
CT).  The  error  was  a  result  of  a  misprint  in  the  latitude  and  longitude  file  of  Weston 
Observatory's  event  location  program  which  resulted  in  an  inaccurate  geographical  placement  of 
station  BCT  on  the  map.  The  erroneous  source  to  receiver  dtetance  value  yielded  group 
velocities  that  were  significantly  higher  than  those  from  the  other  two  regions  and  resulted  in 
what  became  known  as  the  Waterbury  Dispersion  Region  (WDR),  a  region  of  anomalousty  high 
group  velocity. 

After  the  discovery  of  this  error,  Kafka  &  Bowers  (1991)  reprocessed  all  data  from 
paths  that  crossed  the  WDR  using  the  correct  source  to  receiver  distance.  Results  of  this 
correction  indicate  that  Rg  group  velocities  in  the  WDR  are  not  anomalousiy  high,  but  are  on 
average  consistent  with  the  average  group  velocities  for  the  BADR.  For  fiiis  reason,  the  area 
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Tables 


Authors 

Period(sec) 

G.V.(kni/sec) 

Region 

Press  &  Ewing 
(1954) 

8.0-12.0 

3.05 

N.  America 

Bath  (1954) 

3.0-14.0 

3.07 

Europe/Asia 

McEvtIly  & 

Stauder  (1965) 

0.5-5. 0 

1. 5-3.0 

lilinois  Basin/Ozark  Uplift 

Anderson  & 

Dorman  (1973) 

0.2-1. 5 

1. 2-3.1 

Southern  NY/Northem  NJ 

Bath  (1975) 

0.5-1. 8 

2.5-3.0 

Sweden 

Kafka  &  Doliin 

(1985) 

0. 5-2.0 

2. 0-3. 3 

Southern  New  England 

McTigue  (1986) 

0.5-2. 0 

2. 1-2.9 

MAandRI 

Kafka  &  Reiter 

(1987) 

0.4-1. 6 

2.4-3.3 

Southeastern  Mahe 

Tu  (1989) 

0.3-2. 0 

2.3-2.9 

NH  and  Vennont 

this  study 

0.2-2. 2 

2.0-3.29 

MA  and  Southern  NH 
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that  was  designated  as  the  WDR  has  sbioe  been  merged  with  the  quite  extensive  BADR,  and  it  is 
probabiy  rx)  ionger  appropriate  to  use  either  name  to  refer  to  this  dispersion  region.  (Note: 
Figure  1  reflects  the  changes  made  post  Kafka  &  Bowers  (1991).  However .  the  results  of  this 
study  show  that  lateral  variation  can  be  observed  in  the  shallow  crust  beneath  SNE  and  southern 
NH  when  a  denser  station  spacing  is  used. 

5.2  Results  of  This  Study 

The  variation  of  group  velocities  observed  in  this  study  suggests  that  there  is  indeed 
systematic  lateral  variation  in  the  structure  of  the  upper  crust  underlying  SNE.  Various  ways 
of  displaying  the  results  shows  that:  1)  lateral  variations  in  group  velocities  within  the  study 
area  appear  to  be  a  function  of  distance  from  the  San-Vel  quarry,  yet  are  so  subtle  that  they  can 
only  be  observed  on  a  small  scale  and  are  most  easily  observed  when  the  data  are  plotted  in  three 
dimensions.  The  denser  station  spacing  reveals  more  detail  than  what  has  previously  been 
observed;  2)  group  velocity  variations  do  rtot  appear  to  be  a  function  of  azimuth,  therefore, 
there  does  not  seem  to  be  any  evidence  of  lateral  anisotropy  in  the  shallow  cmst  underlying  the 
study  area  or  any  paths  that  coukJ  be  considered  anomalously  fast  or  slow  (excluding  path  NWIO 
which  exhbits  very  high  group  velocities,  however  this  cannot  be  readily  explained  by  the 
geology  or  the  geophysics). 

Figure  15  is  a  graph  of  the  group  velocity  curves  for  all  of  the  paths  considered  in  this 
study.  What  is  most  striking  about  this  graph  is  that  the  majority  of  the  data  falls  within  such  a 
small  range  of  group  velocities.  However,  one  path  (NWIO)  exhfoits  higher  group  velocities. 
This  is  strange  because  the  path  NWIO  does  not  cross  particular  geologic  structures  nor 
terminate  in  any  particular  rock  type  that  would  make  it  distinct  from  the  paths  that  surround 
it,  especially  NW08  and  NW1 1  which  happen  to  exhibit  group  velocities  that  are  typical  of  the 
rest  of  the  study  area. 
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A  few  of  the  group  velocity  curves  appear  to  be  anomalously  slow.  These  paths.  NW03, 
WFM,  and  UXB  doni  appear  to  have  anything  in  common  with  each  other  that  is  obvious  enough 
to  suggest  any  reason  why  they  exhibit  slower  group  velocities  than  the  rest  of  the  group. 

Figure  16  shows  the  average  group  velocity  cunre  for  the  BAOR  from  Kaflta  (1988),  as 
well  as  the  average  and  standard  deviation  for  the  BADR  Including  all  post  1987  data.  The  most 
interesting  thing  about  this  picture  is  the  similarity  between  the  two  curves.  Figure  17  is  a 
graph  that  shows  only  the  group  velocity  cunres  for  data  that  fell  into  categories  1  and  2  (i.e. 
the  best  data).  This  data  was  graphed  to  show  that  even  by  weeding  out  the  poor  quality  data  we 
still  see  the  same  grouping  of  the  cun/es  as  we  did  when  we  Included  all  of  the  data  (Figure  15). 
This  would  suggest  to  me  that  data  quality  has  only  a  limited  effect  on  the  actual  data 
interpretation.  In  other  words  we  can.  in  most  cases,  glean  the  same  information  from  a 
seismogram  of  fair  recording  quality  as  one  of  superior  recording  quality. 

in  order  to  test  for  anisotropy  .plots  were  made  of  group  velocity  vs.  azimuth  for  periods 
between  .5  sec  and  1.2  sec.  Figure  18  shows  the  results.  The  (fiagrams  in  the  left-hand  column 
include  all  of  the  data  used  in  the  study.  The  right-hand  column  represents  only  the  best  quality 
data  (quality  1  and  2  seismograms).  A  third  order  regression  line  was  fitted  to  the  data.  It  is 
apparent  in  these  graphs  that  there  is  very  little  change  in  the  data  amongst  the  different 
periods  and/or  between  the  different  quality  data.  The  variation  of  the  regression  curves  is  very 
slight  and  is  mainly  driven  by  the  three  or  four  anomalous  paths  (NW10,  NW03,  UXB,  and 
WFM)  that  I  discussed  above.  If  I  were  to  remove  these  few  paths,  the  regression  line  would  be 
virtually  flat.  In  either  case,  the  plots  in  this  figure  demonstrate  that  differences  in  azimuth  do 
not  appear  to  have  a  significant  effect  on  group  velocity.  Therefore  we  do  not  observe  any 
evidence  of  lateral  anisotropy  in  the  shallow  crust  beneath  the  study  area. 

Figure  19  shows  the  group  velocity  results  as  a  function  of  distance.  As  in  Figure  18, 
this  figure  displays  each  period  separately  and  also  divides  the  data  into  two  groups.  The  entire 
data  set  is  plotted  in  the  left-hand  column  and  quality  1  and  2  data  are  shown  in  the  right- 
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GROUP  VBJOCiTY 


GROUP  VELOCITY 


Figure  1 6 :  Group  Velocity  averages  and  standard  deviations  for  pre- 
and  post- 1987  data  sets,  inset:  Origin  time  errors  calculated  from 
Figure  14. 
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Figure  1 7 :  Group  velocity  curves  of  quality  1  and  quality  2 
seismograms. 
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Group  velocity  plotted  as  a  function  of  azimuth. 
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Figure  1 8  continued. 
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Figure  18  continued. 
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Figure  1 9  :  Group  velocity  plotted  as  a  function  of  distance. 
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Figure  19  continued. 
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handed  column.  Again,  a  third  order  regression  line  was  fitted  to  the  data,  and  as  in  the  case  of 
azimuth,  the  regression  line  is  virtually  flat.  The  data  show  no  obvious  changes  at  different 
periods  or  between  the  different  quality  data.  Plotted  in  this  manner,  it  appears  that  the  group 
velocity  does  not  strongly  depend  on  distsmce.  However,  there  is  some  slight  incflcation  of 
velocity  increasing  with  distance,  and  when  the  same  data  are  plotted  in  map  view,  we  can  see 
that  the  group  velocities  in  the  area  surrounding  the  San-Vel  quarry  are  noticeably  slower  than 
those  farther  away.  This  strongly  suggests  that  there  actually  are  small  scale  variations  in 
group  velocity.  Figures  20a4),  and  c  are  map  view  plots  of  the  group  velocity  for  three 
different  periods.  The  group  velocities  are  divided  into  four  ranges  and  each  range  is 
represented  by  a  different  symbol. 

Figures  21a,b,  and  c  are  three-dimensional  plots  of  the  same  data  shown  in 
Figures  20a,  b  and  c.  These  figures  also  include  a  contour  map  of  the  same  data  (with  the 
station  array  superimposed)  for  the  same  periods.  These  figures  were  plotted  using  the 
*Surfer”  three-dimensional  plotting  program  published  by  Golden  Software  Inc.,  Golden 
CO.  "Grid",  a  module  within  Surfer,  creates  a  regularly  spaced  grid  from  irregularly 
spaced  data.  An  inverse  distance  inteqx)lation  method  was  used  to  estimate  the  values  of 
group  velocity  for  the  grid  points.  The  output  from  the  Grid  module  was  used  to  create 
the  three  dimensional  and  contour  plots. 

What  is  most  striking  about  these  plots  is  that  with  the  data  plotted  in  three 
dimensions  you  can  see  subtle  undulations  in  the  group  veiocHies  that  are  less  than 
obvious  when  using  other  ways  of  displaying  the  data  (i.e.  Figures  18-  20).  Points  of 
interest  in  Figure  21  are:  1)  the  obvious  spike  at  station  NW10.  Again,  there  doesn't 
seem  to  be  an  obvious  explanation  for  the  anomalously  high  group  velocities  here.  2) 

The  quarry  seems  to  lie  in  a  "trough"  (i.e.  an  area  of  particularty  low  group  velocities). 

With  increasing  distance  from  the  quarry,  group  velocities  systematically  increase 
toward  the  east  and  towards  the  west-northwest,  creating  a  "U-shaped  valley"  of  group 
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velocities  where  the  trough  of  the  valley  appears  to  strike  in  a  northeast  direction.  This 
is  particularly  interesting  because  the  trend  of  this  trough  correlates  well  with  the 
Ciinton-Newbury  Fault  Zone.  One  possible  explanation  of  this  result  is  that  group 
velocities  within  a  fractured  and  faulted  srea  are  likely  to  be  slower  than  group 
velocities  for  Rg  waves  that  travel  through  solid  (unfractured)  bedrock. 
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Figure  20  (a):  Group  velocities  for  Rg  waves  with  a  period  of  0.5 
seconds.  Top  figure  includes  ail  data.  Bottom  figure  includes  only 
seismograms  with  known  origin  times. 
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Figure  20  (b):  Group  velocities  for  Rg  waves  with  a  period  of  0.7 
seconds.  Top  figure  includes  all  data.  Bottom  figure  includes  only 
seismograms  with  known  origin  times. 
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Figure  20  (c):  Group  velocities  for  Rg  waves  with  a  period  of  0.9 
seconds.  Top  figure  includes  all  data.  Bottom  figure  includes  only 
seismograms  with  known  origin  times. 


Figure  21  (a):  Three  dimensional  plot  and  contour  map  of  group 
velocity  for  Rg  waves  with  a  period  of  0.5  seconds  for  the  study 
area.  The  station  array  (dots)  and  the  Ciinton-Newbury  fault  (solid 
line)  are  superimposed  on  the  contour  map 
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Figure  21  (c);  Three  dimensional  plot  and  contour  map  of  group 
velocity  for  Rg  waves  with  a  period  of  0.9  seconds  for  the  study 
area.  The  station  array  (dots)  and  the  Clinton-Newbury  fault  (solid 
line)  are  superimposed  on  the  contour  map 
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Figure  22  is  a  map  of  the  study  area  showing  the  station  array  superimposed  on  the 
regional  geology,  it  is  significant  to  note  that  the  general  trend  of  the  regional  geologic 
structure  as  well  as  the  trend  of  the  Ciinton-Newbury  Fault  Zone  roughiy  parallels  the  trend  of 
the  lower  group  velocity  zone  that  surrounds  the  San-Vei  quarry.  This  figure  shows  which 
paths  lie  within  this  complex  fault  zone  and  which  paths  travel  across  solid  bedrock.  Comparing 
Figure  22  with  Rgure  20  shows  that  the  area  immediately  surrounding  the  San-Vel  quarry 
appears  to  be  characterized  by  slower  group  velocities,  probably  because  the  bedrock  is  highly 
faulted  and  fractured.  Group  velocities  seem  to  systematically  increase  with  increasing  distance 
from  the  fault  zone,  probably  because  the  Rg  waves  spend  more  time  traveling  through  solid 
bedrock  than  fractured  bedrock.  Theref6re,lt  is  concluded  that  based  on  the  results  of  this 
study,  the  shallow  crust  beneath  SNE  and  southern  NH  does  indeed  exhibit  lateral  variations  in 
group  velocity  and  those  variations  seem  to  correlate  to  some  extent  with  the  complex  geologic 
structures  of  that  area.  A  denser  station  spsxiing  and/or  shorter  paths  has  successfuliy  revealed 
more  detail  in  lateral  variation  than  has  been  observed  from  previous  studies. 

As  long  as  there  is  an  observable  Rg  wave  on  the  seismogram,  the  quality  with  which  a 
seismic  event  is  recorded  seems  to  have  very  little  effect  on  the  accuracy  with  which  group 
velocities  can  be  measured.  This  is  most  easily  seen  in  Figures  18  and  19.  It  seems  that  if  an  Rg 
wave  has  been  recorded,  then  a  poor  s^nal-to-noise  ratio  will  not  affect  the  group  velocity 
cunre  that  is  extracted  from  the  seismogram.  However,  it  comes  as  no  suprise  that  if  Rg  energy 
is  absent  from  the  seismogram  to  begin  with,  then  all  the  filtering  in  the  world  will  not  produce 
a  good  group  velocity  curve. 

Finally,  the  shallow  crust  beneath  SNE  and  southern  NH  does  not  display  significant 
lateral  anisotropy  as  does  the  shallow  crust  in  southeastern  Maine. 
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